1

Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to EPSG:3035 and check your changes.
You need to rename one of the id variables or adjust your join accordingly (AGS = district_id).

2

We want a first descriptive visual of the distribution of Covid-19 deaths in Konstanz and the surrounding districts. Calculate the number of Covid-19 deaths (death_rate) by population (population) and multiply with 100,000.

Filter Konstanz (district_id == "08335"), find the surrounding districts, and plot Konstanz and its surrounding districts.
You can use the dplyr function sf::bind_rows() to combine the two spatial objects, “Konstanz” and “Konstanz Surroundings”.

3

Save your data set of Konstanz and its surrounding districts as an ESRI Shapefile.